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Abstract: The discovery of a 125 GeV Higgs boson at the LHC marked a breakthrough in particle physics. 
The relative lightness of the new particle inspires the consideration of a high luminosity Circular Electron 
Positron Collider (CEPC) as a Higgs Factory to study the Higgs boson in a clean environment. At the CEPC, 
the beamstrahlung might represent one of the most important sources of beam-induced backgrounds that will 
impact the detector. It will introduce additional backgrounds to the CEPC detector through the subsequent 
electron-positron pair production and the hadronic process. Therefore its impacts should be carefully evaluated. 
In this paper, the beamstrahlung-induced backgrounds are first estimated with analytical methods and are 
further evaluated in detail with Monte Carlo simulation. The detector occupancy due to the beamstrahlung 
at the location where the first vertex detector layer may be placed is found to be well below 0.5%. Radiation 
levels characterised as non-ionising energy loss (NIEL) and total ionising dose (TID) are estimated to be 
~ 10 11 1 MeV n eg /cm 2 /yr and ~300 kRad/yr, respectively. 
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1 Introduction 


After the Higgs boson was discovered by the AT¬ 
LAS and CMS experiments [l|, 0 at the CERN Large 
Hadron Collider (LHC), it becomes important to 
measure the Higgs boson properties with high pre¬ 
cision, beyond the ultimate reach of the LHC |s|. 
This goal is generally believed to be only achiev¬ 
able at future lepton colliders (so-called “Higgs Fac¬ 
tory”). The proposed Circular Electron Position Col¬ 
lider (CEPC) will primarily operate at the center-of- 
mass energy range of 240 - 250 GeV (near the ZH 
production threshold) with an instantaneous luminos¬ 
ity of 2xl0 34 cm _2 s _1 Q. To achieve the desired high 
luminosity, both electron and positron beams will 
have to be heavily squeezed to very small bunch sizes 
just before collision. At the interaction point (IP), 
trajectories of the electrons or positrons in one bunch 
are bent by the electromagnetic field of the beam par¬ 
ticles of the opposite charge inside a crossing bunch. 
During this process, one special kind of synchrotron 
radiation, called “beamstrahlung” Q, will be emitted. 


Subsequent processes including the electron-positron 
pair production and hadronic event generation Q|, 
both of which involve the energetic beamstrahlung 
photons, are considered among the most important 
detector backgrounds at the CEPC and should be 
carefully evaluated. 


2 Analytical Estimation of the Beam¬ 
strahlung Backgrounds 


The beamstrahlung is usually characterised by the 
beamstrahlung parameter T : 


T __ 2 hw c 
~~ 3 E 


( 1 ) 


where w c = § 7 3 c/p denotes the critical energy of syn¬ 
chrotron radiation, p the bending radius of the parti¬ 
cle trajectory and E the beam particle energy before 
radiation. The higher the T, the more beamstrahlung 
photons with higher energies will be emitted. As¬ 
suming Gaussian charge distributions for the collid¬ 
ing beams, the average T can be estimated with the 
following formula: 
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Table 1. Machine parameters of the LEP2, CEPC and ILC colliders and their average beamstrahlung param¬ 
eters T av calculated with Eq. 0 


Parameters 

Symbol 

LEP2 

CEPC 

ILC250 

ILC500 

Center of mass energy 

E cm [GeV] 

209 

240 

250 

500 

Bunch population 

N [xlO 10 ] 

58 

37.1 

2 

2 

Horizontal beam size at IP 

a x [nm] 

270000 

73700 

729 

474 

Vertical beam size at IP 

Gy [nm] 

3500 

160 

7.7 

5.9 

Bunch length 

a z [/tm] 

16000 

2260 

300 

300 

Horizontal beta function at IP 

/3 X [mm] 

1500 

800 

13 

11 

Vertical beta function at IP 

/3 y [mm] 

50 

1.2 

0.41 

0.48 

Normalized horizontal emittance at IP 

^e x [mm • mrad] 

9.81 

1594.5 

10 

10 

Normalized vertical emittance at IP 

7 Cj t [mm • mrad] 

0.051 

4.79 

0.035 

0.035 

Luminosity 

L [10 34 cm~ 2 s 1 ] 

0.013 

1.8 

0.75 

1.8 

Beamstrahlung parameter 

r a „ [xio- 4 ] 

0.25 

4.7 

200 

620 


5 Nr 2 "i 

T ps- _ e ' _ (o) 

6 a(a x +<Jy)a z 

where r e is the classical electron radius, 7 the Lorentz 
factor of the beam particles, a the fine structure con¬ 
stant, g x /g v the transverse size of the bunch and a z 
the bunch length. 

With the machine parameters of LEP2 [?j], 
CEPC dj and ILC listed in Table [Q the aver¬ 
age beamstrahlung parameters T av of these colliders 
are calculated with Eq. [2] To achieve the required 
high luminosities, the bunch sizes of the CEPC and 
the ILC are expected to be much smaller than that of 
the LEP2, resulting in non-negligible beamstrahlung 
effects at the two machines. 


2.1 Beamstrahlung Photons 

Assuming Gaussian charge distributions for the 
two colliding bunches, the average number of pho¬ 
tons n 7 emitted by each beam particle in one bunch 
crossing (BX) can be approximated as: 


n 7 



1 

(1 + T ^ 3 ) 1 / 2 


(3) 


where A e is the Compton wave length of the electron 
(ij. With the CEPC beam parameters, there will be 
only ~ 0.22 photons emitted by an individual beam 
particle. However the total number of photons per 
bunch crossing can become significant because of the 
high bunch population of 3.71 x 10 11 . 


2.2 Pair Production 


There are two kinds of processes for the pair cre¬ 
ation in the beam-beam interactions, namely coher¬ 
ent and incoherent pair productions. In the coherent 


pair production, real photons convert into electron- 
positron pairs in the strong macroscopic electromag¬ 
netic held of the other crossing bunch. The total num¬ 
ber of coherent pairs in each bunch crossing amounts 
to: 



with 


5(T) 


(7/128) exp(—16/(3T)) (T<1) 

0.295T- 2 C(lnT-2.488) (T»l) 


Eq. [4] suggests that the coherent pair production 
with small T will be exponentially suppressed [101. 
The number of coherent pairs produced per bunch 
crossing at the CEPC is very close to 0. 

In the incoherent pair production, electron- 
positron pairs are produced by the interaction be¬ 
tween two incoming photons, which could be real 
and/or virtual. There are three incoherent processes: 
the Breit-Wheeler process 77—>-e + e _ , in which both 
photons are real; the Bethe-Heitler process e 7 —» 
ee + e“, in which one photon is real and the other is 
virtual; and the Landau-Lifshitz process ee —>• eee + e“, 
in which both photons are virtual. The approximated 
cross sections of the three processes can be calculated 
with the formulae in Refs. (9. 11. T3 . There will be 
roughly 44 pairs from the Breit-Wheeler process, 327 
pairs from the Bethe-Heitler process and 1322 pairs 
from the Landau-Lifshitz process produced in every 
bunch crossing at the CEPC. 

2.3 Hadronic Backgrounds 


In addition to the electron-positron pairs, two col¬ 
liding photons can also produce hadrons with certain 
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Fig. 1. A possible layout of the CEPC interaction region. 


probability. This process is non-negligible because a 
small fraction of the events may contain final state 
particles of high transverse momenta. The cross sec- 
tion of the hadronic process can be parameterised 
as 


13J: 


a H = 
211nb • 


(ah’)' 


+ 297nb 


(ah’)' 


(5) 


where s is the square of the center-of-mass energy 
of two photons. The number of hadronic events per 
bunch crossing at the CEPC will be much less than 

1 . 

The analytical estimation suggests that there will 
be significant amount of beamstrahlung photons per 
bunch crossing at the CEPC. The consequent inco¬ 
herent electron-positron pair production might be the 
dominant background induced by the beamstrahlung, 
while the coherent pair production and hadronic 
events could be negligible. Nevertheless, to fully as¬ 
sess their impacts on the CEPC detector, studies 
based on full Monte Carlo simulation should be per¬ 
formed. 


3 Simulation Studies of the Beam¬ 
strahlung Backgrounds 

The beam-beam interactions have been simulated 
with Guinea-Pig++ (Generator of Unwanted Interac¬ 
tions for Numerical Experiment Analysis Program In¬ 
terfaced to GEANT) [1J], which allows detailed stud¬ 
ies of the emission of the beamstrahlung photons, the 
incoherent pair production and the hadronic events. 
In the simulation, the charge distributions in both 
collision bunches are assumed to be Gaussian. The 


background events generated with Guinea-Pig-t—I- are 
interfaced to Mokka [ll]/GEANT4 0 for full detec¬ 
tor simulation. Hit information such as hit density 
and total ionizing dose are extracted accordingly. 

The primary background particles hitting the 
beampipe or any detector component can cause par¬ 
ticle shower. These secondary particles may enter 
the detector and become important background for 
the detector. Therefore the impact of beam back¬ 
grounds on the detector is highly depending on the 
layout and material budget of the beampipe and the 
detector. Fig. [L] shows the interaction region layout 
used in this study, which features a short focal length 
of L* = 1.5 m, i.e. the distance between the final 
quadrupole QDO and the IP. 

4 Simulation Results 



e [rad] 

Fig. 2. pt versus 8 distribution of the beam¬ 
strahlung photons at the CEPC. 


4.1 Beamstrahlung Photons 

Fig. [2] shows the transverse momentum distribu¬ 
tion versus the polar angle of the beamstrahlung pho¬ 
tons originating from the beam-beam interactions. 
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The photons are confined within a small angle of 
\6\ < 1 mrad, which indicates that most of them 
will leave the interaction region without interacting 
with the beam pipe or any other detector component. 
Therefore they are not the major source of detector 
backgrounds. 

4.2 Incoherent Electron-Positron Pairs 

Unlike the beamstrahlung photons, the electrons 
and positrons from the incoherent pair production 
are usually produced with large polar angle and high 
transverse momentum. They can contribute to the 
detector backgrounds in direct and indirect ways: 

• If the primary particles are with large enough 
polar angle and transverse momentum, they 
will hit the detector directly; 

• If the primary particles are with small polar 
angle, they might still hit the beampipe or 
the downstream detector components and pro¬ 
duce secondary particles, which might be back- 
scattered into the detector. 

Fig. [3] shows the p T versus 9 distribution of the 
electrons and positrons. There is an empty region 
at the bottom-right corner of the distribution plot 
because electrons and positrons with energy below 
5 MeV, which will be constrained in the beam pipe 
by the solenoid magnetic field, are not tracked in the 
simulation with Guinea-Pig-1—h It should be noted 
that the beam pipe and any detector component must 
be kept away from the particles populated in the red 
region to avoid particle showering. It requires addi¬ 
tional optimisations to the beam pipe design and the 
downstream detector components layout in the inter¬ 
action region. 



10 4 10' 3 10 2 10' 1 1 


6 [rad] 

Fig. 3. The pr versus 9 distribution of the elec¬ 
trons and positrons from the incoherent pair 
production. 


In the p T versus p z distribution of the electrons 
and positrons, the major fraction of the particles are 
concentrated in the area under the envelope (so-called 
“kinematic edge”). The developed envelop can be fit¬ 
ted to an imperial formula as the following: 

Pt = 0.0202 p° z 297 (6) 

Assuming perfect helical trajectories in a solenoid 
magnetic filed of 3.5 T, the particles on the envelop 
can develop a trajectory profile as shown in Fig. Q] 
The beam pipe, labeled as the red line, must be kept 
away from the helical trajectories. 



Fig. 4. The helix trajectories of elec¬ 
trons/positrons from the pair production 
in a solenoid magnetic field of 3.5 T. The red 
line indicates the position of the beam pipe, 
which must be kept away from the profile of 
the helical trajectories. 


4.3 Hit Density 

Although most of the particles originating from 
the beamstrahlung will leave the interaction region 
without hitting the beampipe or any detector compo¬ 
nent, a small fraction of them will enter the detector 
directly, or hit the beampipe and/or detector com¬ 
ponents in the forward region and introduce back- 
scattered particles. These primary and secondary 
particles will increase the detector occupancy and 
introduce radiation damage to the silicon detectors 
close to the interaction point. 
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R [mm] 

Fig. 5. Hit density at the VTX and SIT detec¬ 
tor layers 


inc. pairs, p T >5MeV: NIEL, 1-MeV-neutron equlv.flux/cm 2 /yr 



Z [mm] 


Fig. 6. NIEL distribution in the CEPC vertex 
detector. 


The estimated hit densities at the barrel layers 
of the vertex detector (VTX) and the silicon inner 
tracker (SIT) are shown in Fig. [5] On the inner most 
VTX layer at r = 1.6 cm, the hit density amounts to 
0.2 hits/cm 2 /BX. With the design bunch spacing of 
~ 3.6 fj,s for the CEPC, the corresponding detector 
occupancy will be well below 0.5% for the first VTX 
layer with an assumed pixel pitch size of 20 /im and 
a readout time of 20 ps. 

4.4 Radiation Levels 


Although the estimated hit density is low per 
bunch crossing, the detector radiation levels can be 
still considerable given the high repetition frequency 
(~ 2.8 x10 s bunch crossings per second) of the CEPC 
machine. Their impacts on the CEPC detector, in 
particular the vertex detector close to the interaction 
point, should be carefully evaluated. The radiation 
damage in silicon detector can be roughly character¬ 
ized as the non-ionizing energy loss(NIEL) and the 
total ionizing dose(TID). 

NIEL can lead to crystal defects by displacing the 
silicon atom out of their lattice sites (so-called “bulk 


damage”) [17|. The effects induced by any particle 


with a given energy can be normalised to the equiva¬ 
lent damage caused by 1 MeV neutrons. In this study, 
the electron flux at a given point can be obtained by 
tracking all particles with perfect helixes. The NIEL 
per year is calculated with the same method used in 
Ref. fll . assuming machine operation for 10 7 sec¬ 
onds (Snowmass year) and with a safety factor of 10 
taken into account. Fig. [0] shows the NIEL distribu¬ 
tion in the CEPC vertex detector. At the first vertex 
detector layer, the annual value for NIEL is about 
10 11 1 MeV n e9 /cm 2 /yr. 


On the other hand, TID can create ionization at 
the Si-Si0 2 interface and degrade the performance of 
silicon devices 191. TID can be determined with: 


D = 


E, 


dep 


E, 


dep 


P-V 


(7) 


where E dep is the energy deposited in the matter, and 
to, p and V are the mass, the density and the volume 
of the matter, respectively. The annual TID in the 
VTX with full simulation is shown in Fig. [7] and a 
safety factor of 10 is included in the calculation. The 
TID is estimated to be ~300 kRad/yr at the first 
vertex detector layer. 



Fig. 7. TID at the vertex detector layers. 


5 Conclusion 


The backgrounds induced by the beamstrahlung 
have been studied with both analytical and simula¬ 
tion methods. The hit density on the first vertex 
detector layer is about 0.2 hits/cm 2 /BX and the de¬ 
tector occupancy is well below 0.5%. The annual 
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values for NIEL and TID, representing the levels of 
the radiation damage in silicon detector, are esti¬ 
mated to be 10 11 1 MeV n eq / cm 2 and ~300 kRad, 
respectively. Further studies on the beamstrahlung 
effects and other sources of detector backgrounds, as 
well as detailed evaluation of their impacts on physics 
measurements will be pursued in future. 
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